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A  COMPACT  BULK  ACOUSTO-OPTIC  TIME  INTEGRATING  CORRELATOR 

D.A.B.  Fogg 


SUMMARY 

This  report  describes  the  design  and  testing  of  a 
prototype  correlator  using  lead  molybdate  Bragg  cells.  It 
indicates  the  suitability  of  this  design  for  cells  made  of 
materials  with  longer  transition  times,  such  as  tellurium 
dioxide.  Compactness  and  simplicity  are  achieved  by  the 
use  of  a  specially  designed  prism:  An  analysis  of  the  AO, 
TIC  process  and  performance  parameters  is  given. 
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1 .  INTRODUCTION 

Acousto-optical  (AO)  signal  correlators  have  been  described  and  discussed 
extensively  in  past  literature(ref . 1  to  5).  Compactness  is  often  of  prime 
importance  in  such  devices.  This  report  describes  a  dual  beam  bulk  acousto¬ 
optic  time  integrating  correlator  (TIC)  in  which  compactness  and  simplicity 
are  achieved  by  use  of  a  specially  designed  prism.  The  shape  and  transmission 
coating  permit  the  direction  of  laser  beams  of  the  correct  intensity,  angular 
separation  and  polarisation  through  two  Bragg  cells.  Unknown  and  reference 
signals  supplied  to  these  cells  modulate  the  incoming  light.  A  CCD  linear 
photodiode  array  detector  is  used  to  sum  and  square  the  modulated  light.  The 
detector  output  contains  the  desired  correlation. 

Correlator  analysis  is  usually  simplified  by  using  analytic  signal  equivalents 
and  omitting  the  optical  beam  expression.  A  fuller  analysis  is  given  in 
Section  2.  Although  more  laborious,  it  leads  to  a  better  understanding  of  the 
process . 

Specifications  for  the  prism  vary  with  Bragg  cell  type.  Prism  parameters  are 
described  analytically,  enabling  practical  constructuion  details  to  be 
calculated  to  suit  any  Bragg  cell. 

The  correlator  architecture  is  described.  A  prototype  was  made  using 
available  lead  molybdate  cells  and  experimental  results  are  given  for  the 
correlation  of  a  number  of  signals.  Performance  criteria  are  discussed  in 
Section  6,  and  the  bases  for  these  criteria  are  reviewed  in  Appendix  II. 

2.  TIC  ARCHITECTURE 

In  bulk  AO  TIC  correlators,  the  deflected  first  order  beam  of  one  cell  may  be 
spread  and  introducted  at  the  Bragg  angle  to  the  second  cell.  The  deflected 
beam  from  the  second  cell  contains  the  information  which  on  integration  at  the 
detector,  provides  the  desired  correlation.  Alternatively,  separate  first 
order  beams  from  each  cell  can  be  produced  and  combined  by  reflection  or 
refraction.  Surface  acoustic  wave  devices  have  been  designed(ref . 7)  to 
produce  two  collinear  first  order  beams.  The  architecture  developed  here 
produces  collinear  first  order  beams  from  bulk  AO  devices  using  a  specially 
designed  prism,  resulting  in  an  efficient  compact  assembly. 

If  the  direction  of  the  incident  laser  beam  to  a  Bragg  cell  has  a  component  in 
the  direction  of  propagation  of  the  acoustic  strain  wave  in  the  cell,  the 
deflected  light  is  downshifted  in  frequency  and  is  referred  to  as  the 
-1th  order(ref . 3) .  Alternatively,  if  the  incident  beam  has  a  component  in  the 
opposite  direction  to  the  acoustic  strain  wave  a  +lth  order  (upshifted) 
deflected  beam  is  obtained.  By  orientating  two  Bragg  cells,  and  choosing 
appropriate  acoustic  propagation  directions  it  is  possible  to  obtain 
identically  shifted  collinear  deflected  beams,  one  from  each  cell.  Figure  1 
shows  the  required  arrangement  to  achieve  either  upshifted  or  downshifted 
collinear  first  orders.  In  both  cases,  the  angle  between  the  cell  axes  is 
twice  the  Bragg  angle  and  the  two  incident  beams  have  a  divergence  of  four 
times  the  Bragg  angle. 

The  deflected  +1^  and  the  undeflected  0^  beams  from  the  first  cell  pass 

through  the  second  cell,  B.  The  0^  beam  makes  an  angle  30g  with  cell  B.  This 

is  not  close  enough  to  the  Bragg  angle  to  give  rise  to  further  deflected  beams 
and  it  is  transmitted  virtually  unchanged.  The  +1^  beam  incident  on  cell  B  at 

0g  gives  rise  to  a  deflected  beam  coincident  with  0g.  It  suffers  a  small 

intensity  loss  equal  to  the  deflection  efficiency.  (Of  the  order  of  6%  for 
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the  lead  molybdate  cells  and  the  RF  power  used.) 

The  correlator  is  shown  schematically  in  figure  2.  A  beam  from  a  10  mW  laser 
was  spatially  filtered  and  expanded  to  provide  a  10  mm  wide  by  1  mm  high  beam. 
P-polarisation  was  used  (ie  the  electric  vector  in  the  plane  of  incidence)  to 
optimise  the  performance  of  the  Bragg  cells.  A  phase  equalising  plate  was 
used  in  one  half  of  the  beam.  This  plate  could  be  rotated  5°  to  obtain  a  full 
360°  phase  change.  The  prism  produced  two  equal  intensity  10  mm  x  1  mm  beams 
diverging  at  an  angle  of  48  but  precisely  in  the  same  horizontal  plane.  The 

D 

undeflected  beams  were  stopped  shortly  after  separating  from  the  coincident 
first  orders.  A  1024  element  CCD  array  was  used  as  the  det'.ctor. 


3.  TIC  ANALYSIS 

The  RF  signal  modulated  in  amplitude  and  phase  and  directed  to  the  Bragg  cell 
transducer  can  be  denoted  by, 


s(t)  =  A(t)  cos{2tt  f  t  +  a(t)} 

c 


(1) 


=  Re  {A(t)  exp  j  (2irfct  +  a(t)]} 


(2) 


where  Re  denotes  "the  real  part  of". 

Within  the  cell  this  driving  waveform  is  given  by(ref.3)  s(t  -  x/v)  or 
s(t  +  x/v  -T)  where  v  is  the  acoustic  velocity  within  the  cell  and  T  is  the 
transit  time  (cell  width  divided  by  acoustic  velocity).  The  former  expression 
represents  sound  wave  propagation  in  the  positive  direction  and  the  latter  the 
reversed  propagation. 

The  object  is  to  find  the  autocorrelation  function  ^ss(i)  of  the  signal  s(t). 

It  is  possible  to  simplify  the  analysis  by  considering  the  signal's  associated 
analytic  function  given  by 


z(t)  =  s(t)  +  j  s(t) 


or 


z(t)  =  A(t)  exp  j  [  2ir  f  ^  t  +  a(t)] 


(3) 


for  s(t)  given  by  equation  (2) 

Here  s(t)  is  the  Hilbert  transform  of  s(t)(ref.8).  The  envelope  e(t)  of  z(t) 
is  obtained  by  multiplying  by  exp  (-j2irf  t)  (ie  extracting  the  carrier)  so 

that 


E(t)  =  z(t)exp(- j2nf ct) 


(4) 


The  Fourier  transforms  S(w),  Z(w)  and  E(w)  of  s(t),  z(t)  and  e(t)  respectively 
are  related  by(ref.8) 
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Z(w) 


|  2  S(w)  w  >  0 


0  w  <  0 


(5) 


and 


E(w)  =  Z  (w  +  w  ) ,  u  =2irf 
o  o  c 

Therefore  if  the  spectrum  E(w)  of  the  envelope  is  confined  to  the  positive 
region,  e(t)  contains  the  information  concerning  the  slow  amplitude  and  phase 
variations  of  the  signal  s(t)  and  satisfies  the  intuitive  concept  of  an 
envelope  (ie  the  carrier  frequency  f  is  greater  than  the  highest  frequencies 

occuring  in  the  information  component  of  the  signal). 

It  is  not  sufficient  to  find  the  autocorrelation  function  R  (t)  of  this 

ee 

envelope  since  it  is  not  equal  to  the  autocorrelation  function  Rss(t)  °f  the 
real  signal.  They  are  related  by 


Ree(t)  =  Rzz(t)exp(-j  2irfc  t)  =  2  [Rgs  (i )  +  j  R^  (t )  ]  (6) 


where  R  (t)  is  the  Hilbert  transform  of  R  (t)(ref.  8).  Therefore 
ss  ss 


Rss(t)  =  i  Re  {Re£(i)}  (7) 

A  number  of  authors (ref . 3)  for  example  treat  the  case  where  the  analytic 
signal  z(t)  or  a  complex  equivalent  to  s(t)  is  assumed  to  be  the  input  to  the 
TIC.  The  output  is  usually  shown  as, 

/•t  i 

0 (t)  =  I  lijz(t  -  x/v)  +  £  jz(t  +  x/v  -T) | 2dt  (8) 

•  t 

o 


Expanding,  this  becomes 

0(t)  = 

il 

rti 

z*(t  -  x/v)dt  +  z2(t  +  x/v  -T)dt 

( 

J 

t 

0 

t 

0 

+ 

* 

•t  i 

Re{z(t  -  x/v)z*(t  +  x/v  -T)}dt 

J 

t 

o 

where  t|  ■  t  is  the  integration  time. 

The  last  term  is  £  Re{R  (2x/v  -T)}.  From  equation  (6) 
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(exp  j2irf  i)  R  (2x/v  -T)  =  2[R  (2x/v  -T)  +  j  R  (2x/v  -T)  ]  (9) 

C  22  S  5 


so  that 


i  Re  {R2Z(2x/v  -T)>  =  Rgs (2x/v  -T)  (10) 

It  is  then  implied  that  using  the  real  signal  instead  of  the  analytic  one  will 
produce  the  same  term. 

It  is  instructive  to  consider  this  process  in  more  detail,  incorporating  the 
direction  of  propagation  of  the  light  and  delaying  use  of  analytical  functions 
until  after  an  expression  for  the  output  light  is  obtained.  The  laser  beam 
incident  on  one  face  of  a  Bragg  cell  can  be  represented  by  the  coherent  plane 
wave 


L(£,d^t)  =  Re  {a  exp  j  2ir(vt  -  1/ X  jr.d)}  (11) 

where  r  =  xi  +  zk  is  the  position  vector  of  a  point  P  in  space  (see  figure  3) 

d  =  sin  0i  +  cos  8  k  is  a  unit  vector  in  a  fixed  direction. 

and  X  is  the  wavelength  of  light  in  the  medium  of  propagation.  The 
quantity  L(r,d,t)  represents  a  plane  wave  since  it  is  constant  over  the 
plane  r.d  at  any  time  t. 

The  optical  output  of  the  Bragg  cell  is  determined  by  assigning  to  it  a 
transmittance  tr(x,t)  equivalent  to  that  of  a  sinusoidal  phase  grating  (see 
reference  9) 

tr(x,t)  =  (Cq  exp  js(t  -  x/v)  (12) 

where  Cq  is  a  function  of  the  incidence  angle  0^  and  represents  the  strength 
of  interaction.  The  defelcted  output  U(x,t)  is  therefore  using  equation  (2), 


U(x,t)  =  Re  {a  exp  j  2ir(vt  -  1 / X  r.d'))  Cq  exp  js(t  -  x/v) 

where  d'  is  d  with  the  sign  of  the  ith  component  chosen  to  account  for  the 
deflection  of  the  light  beam.  Expanding  the  second  exponential  and  taking  the 
first  two  terms  this  becomes 


U(x,t)  =  Re  {a  exp  j  2ir(vt  -  1  / X  r.d'))  Cq { 1  +  js(t  -  x/v)J 


The  deflected  light  U'(x,t)  is  represented  by  the  second  term  in  the  square 
brackets.  This  can  be  written 


=  Re  {a  exp  j  2u(vt  -  1/X  r.d')} 


j  Re  _{A(t_j- jt/v)  exp  j{2xfe(t  -  x/v)  +  a(t  -  x/v)]) 


U'(x.t) 
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Substituting  for  £  and  d  and  using  the  relationship 


Re  {exp  j  0  x  >  Re  {exp  j02}  =  i  Re  {exp  j(81  +  82)  +  exp  j  C  8  x  -  02)} 


this  becomes , 


U'(x,t)  =  i  aA(t  -  x/v)  Re  'C  exp  j 


2ir(v  +  fc)t  +  a(t  -  x/v) 


2ir 


/ 

'nx 


i  i  sin  0.  + -0^  cos  8.  -  f  x/v 

\A  1  A  1C 

o  o 


+  CQ  exp  j 


2ir(v  -  fc)t  ■  a(t  -  x/v) 


+  2u(-^-  sin  0  .  -  cos  0.  +  f  x/v 
\A  1  Xn  i  c 

o  ° 


Here  A  =  nX ,  X  is  the  light  wavelength  in  vacuo  and  n  is  the  refractive 
o  o 


index  of  the  Bragg  cell  material.  This  represents  two  sidebands,  one  at  an 
upshifted  optical  frequency  (v+f)  and  the  other  downshifted  to  (v-f).  When 


the  angle  of  incidence  0.  is  ±0_,  the  Bragg  angle,  one  of  these  sidebands  is 

1  D 

maximised  and  the  other  disappears  (Cq  =  0  or  =  0). 


Consider  the  upshifted  sideband  U2(x,t)  is  equation  (5).  This  has  the  form 


Ul(x,t)  =  i  j  aA(t  -  x/v)  cos[2tt(v  +  f  )t  +  a(t  -  x/v)  +  C^x)]  (13) 


where 


C  i  (x) 


nx 

o  i  n  z 

0  .  +  T—  cos 

V*  sin 

X 

i  X 

o 

o 

0  . 
1 


Notice  that  this  is  not  identical  to  s(t  -  x/v).  The  carrier  is  now  v  +  f 

c 

and  there  is  an  additional  phase  term  Cx(x)  related  to  the  spatial  progression 
of  the  signal  across  the  cell. 

From  the  other  cell  the  output  can  be  shown  to  be 


U2(x,t)  =  i  j  aA(t  +  x/v  -T)  cos[2tt(v  +  f£)t  +  o(t  +  x/v  -T)  +  C2(x)  (14) 

where 


C2(x) 


2v 


^  sin  0, 
X  : 


-t—  cos  0.  +  f  x/v  -  f  T 
X  ic  c 

o 


The  output  Ui  and  U2  are  added  coherently  and  the  intensity  registered  by  the 
detector  is  the  square  of  this  sum.  This  is  integrated  for  a  time  ti  -  tQ, 
producing  an  output  voltage  proportional  to 
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f1 1 

I  lUjCx.t)  +  U2(x,t) | 2dt 
J  t 


Since  Ux  and  U2  contain  phase  terms  it  is  difficult  to  eliminate  the  carrier 

The  best  way  of  interpreting  equation  (15)  is  via  analytic  functions.  Th 
analytic  function  associated  with  Ux(x,t)  is  given  by 

zx(t  -  x/v)  =  i  aA(t  -  x/v)exp  j  [ 2'tt ( v  +  f^)  +  a(t  -  x/v)  +  Cx(x)]  (16 

A  corresponding  expression  exists  for  U2.  Replacing  Ux  and  U2  i 
equation  (15)  by  zx  and  z2 


I  =  J  |Zi(t  -  x/v)  +  z2(t  +  x/v  -T) | 2dt 


but  from  equation  (3) 


z(t  -  x/v)  =  A(t  -  x/v)  exp  j[2it  f^(t  -  x/v)  +  o(t  -  x/v) 


so  that  from  equation  (16) 


Zj(t  -  x/v)  =  1  j  a  exp  j  2irvt  exp  j  Cj(x)z(t  - 


where 


C  i  (x)  = 


C  x ( x )  +  fc  x/v 


and  similarly 


z2(t  +  x/v  -T)  =  i  j  a  exp  j  2irvt  exp  j  C2(x)z(t  +  x/  v  -T) 


where 


C2' (x)  =  C2(x)  -  f  x/v  +  f  T 


Therefore  I  becomes 
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I  -  i  a2  |exp  j  2nvtj2  |z(t  -  x/v)  exp  Ci(x) 

J  t 


+  z(t  +  x/v  -T)  exp  j  C2(x)|2dt 


Since  | exp  j 9 1  =  1  for  real  9,  this  can  be  expanded  as 


ftl 

I  =  i  a2  /  |  z(t  -  x/v) | 2dt  +  i  a2  j  |z(t  +  x/v  -T)|2dt 

Jt  „  t 

o  o 


ftl 

+  £  a  Re  z(t  -  x/v)z*(t  +  x 


-  t 


/v  -T)dt 


(18) 


The  last  integral  has  already  been  considered  following  equation  (9),  and  is 
the  desired  autocorrelation  function  R  g  (2x/v  -T)  multiplied  by  the  light 

intensity  constant,  a.  If  the  integration  time  is  large  compared  with  1  / f c 

the  first  two  integrals  of  equation  (18)  are  constant  depending  only  on  the 
integration  time  ( t x  -  t0).  The  expression  equation  (18)  is  therefore  of  the 

form 


I  =  bias  +  R^^ (2x/v  -T)  (19) 

This  result  is  also  produced  by  the  integral  in  equation  (15).  Thus  the  TIC 
is  seen  to  produce  the  desired  autocorrelation  function  in  a  time  frame 
compressed  by  a  factor  of  2.  Both  amplitude  and  phase  modulation  in  the 
signal  are  seen  to  be  preserved  during  processing  by  the  Bragg  cell.  The 
initial  RF  carrier  of  frequency  f  is  replaced  by  a  light  carrier  of  frequency 

i'  +  f  (or  v  -  f  ).  The  correlation  function  is  then  formed  by  coherent 
c  c 

summation,  squaring  and  integration  in  the  detector.  The  analysis  shows  that 
the  correlator  can  be  used  with  pure  amplitude  modulated,  (a(t)=0),  pure  phase 
modulated  signals  (A(t)=l)  or  both  amplitude  and  phase  modulated  signals. 
This  is  discussed  further  in  Section  5. 


4.  THE  PRISM 

The  prism  is  required  to  superimpose  two  parts  of  a  polarised  spread  laser 
beam  each  10  mm  horizontally  by  1  mm  giving  rise  to  two  beams  with  a  resultant 
divergence  of  40g,  and  equal  intensities.  One  beam  (see  figures  11  and  4) 

enters  surface  3  at  normal  incidence,  is  reflected  from  surface  2  and  emerges 
from  surface  1.  The  second  beam  undergoes  a  reflection  from  surface  1. 

4.1  Design  geometry 

The  prism  is  defined  by  any  four  of  the  angles  given  below,  where  the  angle 
between  surface  i  and  j  is  denoted  by  ajj  (see  Appendix  I) 
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“12 

=  Kn/4  + 

<»U 

=  3  tt/4  - 

On 

=  s in  1  [ : 

a  j  2 

=  |  sin  1 

0)5 

=  tt/2 

-  8  - 

8g)  *  i  sin  1  [sin (it/ 4  -  30g)/n] 
9B 

in(ir/ 4  -  30g)/n) 

[  1/n  sin(TT/4  -  39g)  ]  +  ir/8  +  ±0fi 


(20) 


Here  n  is  the  refractive  index  of  the  glass  used  in  the  prism. 

A  prism  was  manufactured  in  the  DRCS  optical  laboratories  for  use  with  lead 
molybdate  Bragg  cells.  Type  BK7  glass  was  used  and  the  surfaces  were 
ground  and  polished  to  reduce  wavefront  distortion  at  each  surface  to 
better  than  X/10  for  X  =  632.8  nm.  Since  the  Bragg  cell  output  is 
sensitive  to  the  incident  angle  (see  figure  5)  an  attempt  was  made  to  make 
the  prism  angles  as  precisely  as  possible.  The  actual  angles  achieved  were 
measured  using  a  goniometer  and  were  found  to  be  within  2"  of  the  desired 
values.  Pyramidal  error  was  of  prime  concern.  No  pyramidal  error  could  be 
observed  in  the  final  product.  Details  of  the  phase  corrector  plate  used 
to  compensate  for  the  effects  of  path  differences  between  the  two  beams  are 
given  in  Appendix  I . 

4.2  The  effect  of  rotation 

In  order  to  accommodate  variation  in  Bragg  angle  and  slight  errors  in 
achieving  the  required  prism  angles,  the  prism  can  be  rotated.  This 
changes  the  divergence  of  the  two  outgoing  beams.  The  relation  between 
prism  rotation  A0  and  beam  divergence  0  is  derived  in  Appendix  I. 

For  small  A0  this  reduces  approximately  to  a  simple  linear  relationship 


0  =  C  A0  +  Cj 

o 


(21) 


For  lead  molybdate  C  =  0.78  (±0.02)r  and  Cj  =  0.028r  (=40a).  The  prism 

O  D 

can  also  be  placed  on  a  gonimeter  and  rotated  relative  to  the  incident  beam 
until  the  output  beam  divergence  is  precisely  40  . 

D 

A  mathematical  model  of  the  Bragg  cell(ref.ll)  based  on  the  Klein  and  Cook 
analysis(ref . 10)  was  used  to  produce  a  graph  of  the  relationship  between 
output  intensity  and  angle  of  incidence  for  the  lead  molybdate  cell  using  a 
wavelength  of  X  =  632.8  nm.  This  is  shown  in  figure  5  together  with 
experimental  values. 

This  function  is  the  Fourier  transform  of  the  rectangular  aperture  of  the 
incident  beam.  At  the  peak  a  rotation  of  6a  =  1.4  mrad  (5  ft)  produces  a 
detectable  decrease  in  intensity.  To  avoid  reducing  the  sensitivity  of  the 
system  the  prism  manufacturing  tolerances  need  to  ensure  that  the 
divergence  is  within  ±5a  of  20„.  The  actual  error  was  0.02  mrad  (4")  and 

D 

in  the  next  section  it  is  shown  that  this  error  can  be  made  arbitrarily 
small  by  rotation  of  the  prism.  The  prism  therefore  does  not  limit  the 
performance  of  the  system. 
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4.3  Crossover  distance 

To  achieve  the  desired  collinearity  of  the  deflected  beams  the  crossover 
point  of  the  two  input  beams  should  be  between  the  two  Bragg  cells.  The 
crossover  point  can  be  adjusted  by  varying  the  relative  position  of  the 
beams  entering  the  prism.  This  point  can  readily  be  determined  by  trial 
and  error  once  the  prism  has  been  manufactured.  However,  it  is  preferable 
to  ensure  that  the  appropriate  crossover  distance  is  achievable  when 
specifying  the  prism  dimensions  before  manufacture.  The  crossover  distance 
d^  of  two  rays,  one  in  each  beam  is  given  by,  (see  figure  13) 

dc  =  sin(£  -  20b){ (xd  -  xQ)2  +  (YQ  -  YG ) 2 > ^ / ( 2  sin  20fi) 

where 


S  =  it  -  dji,  Z  =  a3  2 

Xp  =  {(Sj  -  Si)tan2f  +  (S3  -  S1)cot2K  +  S3tan£}/ (cot2?f  +  tan£) 
Yd  =  (xD  -  S3)tan£ 

Xq  =  S3  +  S2 

and 


Yq  =  S2tan?  . 

For  a  given  prism  this  can  be  reduced  to  a  simple  relationship  namely, 


d  —  Cj  S2  +  C2 
c 


where  C3  and  C2  are  constants. 

4.4  Coating 

Surfaces  1,  2  and  3  were  coated  to  satisfy  relative  beam  intensity 

requirements.  Surface  4  was  coated  for  maximum  transmission  out  of  the 
prism.  Surfaces  5,  top  and  bottom  were  not  polished. 

Denoting  the  transmittance  and  reflectance  of  the  ith  prism  surface  by  t_^ 

and  r^  respectively,  the  design  requirement  of  equal  beam  intensity  is 

t3r2t3  =  r3,  given  n,  X,  and  p-polar isation .  Surface  3  was  antireflection 

coated,  (t3  =  1).  Because  the  angle  of  incidence  at  surface  2  was  almost 

the  Brewster  angle,  a  dielectric  stack  was  impractical.  An  aluminium 
coating  was  used  and  therefore  r2  =  0 . 83(ref . 12) .  For  a  dielectric  stack 

r  +  t  =  1.  It  follows  that  t3  =  0.55  and  r3  =  0.45.  Coating  details  are 
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given  in  Table  1.  The  measured  beam  intensities  of  the  coated  prism  were 
within  1%. 


TABLE  1.  PRISM  DIELECTRIC  COATINGS 


Surface 

Layer 

Material 

Opt ical 
thickness , 
(nm) 

3 

1 

Zr02 

380 

2 

MgF  2 

140 

1 

1,3,5 

MgF, 

180 

2,4.6 

Zr02 

180 

5.  EXPERIMENTAL  RESULTS 

TIC  correlators  can  process  a  variety  of  signals  including  CW,  linear  FM  and 
direct  sequence  types ( ref . 7 ) .  It  is  instructive  to  use  both  amplitude 
modulated  and  phase  modulated  signals  to  test  the  system.  As  examples  of  pure 
amplitude  modulated  signals,  on-off-keyed  pulsed  CW  (00K/CW)  signals  were 
used.  These  produce  regular  correlation  pulses  across  the  detector  indicating 
the  degree  of  apodisation.  As  examples  of  pure  phase  modulated  signals,  on- 
off-keyed  binary  phase  shifted  (OOK/BPSK)  and  spread  spectrum  binary  phase 
shift  keyed  (SPREAD/BPSK)  signals  were  used.  The  former  for  comparison  with 
the  OOK/CW  signals  and  the  latter  as  examples  of  single  correlation  peak 
output . 

For  these  cases  the  input  signals  s(t)  given  in  general  by  equation  (1) 
be  come (ref. 13) 


OOK/CW  sx  (t) 


rect 


cos  2irf  t 
c 


(22) 


OOK/BPSK  s 2  (t) 


N/2 

V 

L „ 

n=-N/2 


{2  s  (t  -  2nt  )  -  s  (t  -  nt  )) 
o  p  o  p 


(23) 


SPREAD/BPSK  s,  (t) 


N/2 


n=-N/2 


1)  s  (t  -  nt  ) 

°  P 


(24) 
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where  s  (t)  =  rect  t  sin  2irf  t,  f  <<  1/t  ,  -t  /2  <  t  <  t  / 2 ,  and  t  is  the 
o  c  c  p  p  P  p 

pulse  duration.  The  terms  a  are  either  zero  or  one  and  form  a  linear 

n 

recursive  sequence  (LRS).  s^t)  was  synthenised  using  a  signal  generator  and 

a  device  was  built  to  phase  shift  a  suitable  carrier  to  obtain  s2(t).  For  the 

spread  spectrum,  pseudo  noise  signal  s3(t),  a  127  bit  (R7)  LRS  code  was 

generated  using  a  shift  register  with  appropriate  feedback.  The  LRS  generator 

was  defined  by  the  polynomial  x7  +  x1  +  x°  =  0.  By  adding  the  output 

(modulo  2)  from  the  last  two  register  positions  an  identical  LRS  with  lag  7 
was  obtained  and  used  to  generate  the  reference  pseudonoise  signal. 

The  circuit  design  allowed  for  modulation  rates  of  f^/2S,  s  =  1,2,...,  7.  The 

observable  lag  between  the  reference  and  test  signals  was  restricted  to  the 
cell  transit  time  divided  by  the  pulse  duration  t/t  (see  Table  3). 

P 

The  normalised  autocorrelation  functions  for  the  above  signals  are(ref,13): 


OOK/CW  and  BPSK  R  (t) 
s  i  s  j 


R  (i) 
s2s2 


n=-N/2 


t-2nt  /  .  \ 

_ £.1  -  ill)  (25) 

t  '  t  / 

P  P 


III  <  t 


SPREAD/BPSK 


R  (t)  = 

S  3  S  j 


T  >  t 


The  correlation  peaks  are  triangular  in  shape.  The  00K  signals  produced  peaks 
at  every  second  bit  position  across  the  cell,  (see  figure  8).  The  spread 
signals  each  produced  a  simple  correlation  peak  one  bit  wide,  readily 
discernable  above  the  pedestal,  (see  figure  9).  The  pedestal  showed 
indications  of  aperture  diffraction,  the  Gaussian  laser  profile  and  a  periodic 
interference  due  to  a  transparent  cover  over  the  detector  elements.  The 
latter  could  be  removed  by  using  a  detector  without  a  sealing  window. 
Pedestal  subtraction  can  be  used  to  largely  remove  the  effect  of  the  pedestal. 
Figure  9(a)  shows  the  normal  correlator  output  peak.  Figure  9(b)  shows  the 
output  when  a  180°  phase  delay  was  incorporated  into  the  reference  RF  signal. 
Subtracting  the  two  removes  the  pedestal  and  increases  the  height  of  the 
correlation  peak.  Subtraction  was  implemented  by  digitising  the  output  of  the 
1024  element  linear  photodiode  array  using  a  12  bit  analog-to-digital 
converter.  The  conversion  time  was  2.8  ps  giving  a  3  ms  scan  rate.  The 
results  were  stored  on  HP85  computer  disc  files.  The  results  of  the  pedestal 
subtraction  are  shown  in  figure  10.  With  better  beam  shaping  optics  and  a 
windowless  detector  the  base  would  be  as  flat  as  the  outer  edges  of 
figure  10(c). 

The  test  cases  run  and  the  results  observed  are  shown  in  Table  2. 


TABLE  2.  TEST  RESULTS 


Sw i tch 

(s) 

Division 

factor 

(2s) 

Modu l at  ion 
rate. 
fi  =  fc/2s 

(MHz) 

Pul  se 
durat  i  on 
t  =1/2° 
p 

ifiS) 

Ch  i  p 
rate 

T  =1/t 
c  p 

(MHz) 

Max i mum 
lag 

( pu l ses 
in  cell) 

Correlation 

peak* 

BPSK/ 

OOK 

(mV) 

BPSK 

SPREAD 

(mV) 

Pu l sed 
cw 
(mV) 

H 

4 

19.5 

0.026 

105 

>200 

>100 

>400 

H 

8 

9.8 

0.05 

52 

>200 

>400 

>400 

n 

16 

4.9 

0.  10 

26 

>200 

>300 

>300 

BB 

32 

2.4 

0.21 

12 

>200 

>300 

>300 

m 

64 

1.2 

0.41 

6 

>200 

NA 

NA 

■ 

128 

0.61 

0.82 

1.2 

3 

- 

NA 

NA 

f  =  78  MHz,  transit  time  of  cell  =  2.7  /iS 

*  See  figures  8,  9  and  10  to  obtain  an  indication  of  S/N  ratio 


6 .  PERFORMANCE 


The  input  signal  to  a  correlator  will  be  contaminated  by  noise  so  that  the 
output,  equation  (13),  will  be  of  the  form 


c  l 

| U : ( x , t )  +  n(t)  +  U2 (x , t ) | 2dt 

*2 


(27) 


In  a  correlator,  the  energy  from  a  spread  signal  is  concentrated  into  a  narrow 
output  pulse.  The  total  energy  of  the  noise  is  unchanged  by  the  correlator, 
but  spread  in  frequency  since  the  spectrum  of  the  noise  X  signal  product  is 
the  convolution  of  the  noise  and  signal  spectra.  The  noise  in  the  signal 
bandwidth  is  reduced.  Figure  6(a)  shows  the  changes  in  signal  and  noise 
spectra  brought  about  by  the  correlator,  for  the  case  of  a  spread  spectrum 
signal  and  broadband  noise.  Figure  6(b)  shows  the  narrow  band  noise  case. 

Since  AO  correlators  are  essentially  equivalent  to  matched  filters,  the  well 
known  assessments  of  matched  filter  performance  are  appl icable(ref . 14, 15) . 
The  enhancement  of  signal  to  noise  ratio  (ie  processing  gain)  achieved  is 
given  by  the  product  of  the  effective  duration  of  the  signal  correlated  and 
its  bandwidth.  Review  of  classical  and  current  estimates  of  processing  gain 
are  given  in  Appendix  II. 


The  processing  gain  PG  of  a  spatial  integrating  correlator  is  limited  to  the 
product,  Bx ,  of  the  Bragg  cell  bandwidth  B,  and  its  transit  time  x.  In  a  time 
integrating  correlator,  the  factor  limiting  PG  is  either  the  dynamic  range  R 
of  the  detector,  or  the  inherent  PG  of  the  signal,  whichever  is  the  smallest. 
In  the  former  case,  it  can  be  shown(ref.l)  that  the  "effective"  time  bandwidth 
product  for  a  TIC  is  R2 .  In  a  TIC  the  transit  time  limits  the  observable  lag 
between  reference  and  test  signals.  The  maximum  observable  lag  is  the  cell 
transit  time  divided  by  the  pulse  duration,  t  .  For  a  10  MHz  chip  rate 


t  =  0.1  us. 
P 


For  the  prototype  correlator  t  =  3  ps  allowing  a  maximum  lag  of 


30  pulses.  In  many  applications  much  larger  maximum  lags  are  required 
necessitating  Bragg  cells  with  larger  apertures.  Tellurium  dioxide  cells 
currently  offer  the  largest  apertures  at  low  frequencies  when  operated  in  the 
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slow  shear  propagation  mode.  ERL  has  recently  purchased  two  types  of 
tellurium  dioxide  cells  with  apertures  of  50  ps  and  70  ps  and  it  is  proposed 
to  construct  correlators  using  these  cells  which  will  have  maximum  lags  of  500 
and  700  pulses  respectively.  Dynamic  range,  including  the  Bragg  cell's 
inherent  dynamic  range,  is  discussed  in  more  detail  in  Appendix  II. 


TABLE  3.  PERFORMANCE  PARAMETERS 


Bragg  cell 
type 

Frequency  range 
(MHz) 

fmin  cf  fmax 

Bandwidth 

(MHz) 

Transit 

time 

(s) 

Time 

bandwidth 

product 

Resolution 

l_ 

T 

Maximum 
lag  at 

10  MHz 
chip  rate 
(pulses) 

Lead 

Molybdate 

Tellurium 

60 

80 

100 

40 

3 

120 

333  kHz 

30 

dioxide 

deflector 

Tellurium 

30 

45 

60 

30 

50 

1500 

20  kHz 

500 

dioxide 

deflector 

35 

50 

65 

30 

70 

2100 

14  kHz 

700 

Detector  dynamic  range  500:1 


Frequency  resolution  6f  is  determined  by  the  width  D  of  the  light  beam. 
6f  =  v/D  where  v  is  the  velocity  of  sound  in  the  cell.  Resolution  is  best 
when  D  is  the  width  of  the  cell  and  then  6f  =  1/t.  Typical  values  for  the 
lead  molybdate  cell  used  in  the  prototype  TIC,  and  the  proposed  tellurium 
dioxide  correlator  (50  ps  or  70  ps)  are  shown  in  Table  3. 


7.  CONCLUSIONS 

In  this  report,  a  specific  bulk  A0  TIC  design  has  been  described.  It  features 
collinear  deflected  first  order  beams.  This  has  a  number  of  distinct 
advantages.  It  ensures  coplanar  optical  wavefronts  at  the  detector 
eliminating  the  Young's  fringes  seen  in  some  designs.  It  also  avoids  the  need 
for  2  additional  lenses  and  a  frequency  plane  stop  after  the  Bragg  cells  to 
eliminate  the  undeflected  beams.  In  this  architecture,  the  deflected  and 
undeflected  beams  are  separated  by  their  natural  divergent  propagation.  No 
similar  correlator  architecture  has  been  reported  in  the  literature  on  bulk  AO 
time  integrating  correlators. 

The  prism  used  in  the  correlator  is  also  novel  and  facilitates  the  directing 
and  intensity  control  of  the  correlator's  input  laser  beams.  Analytic 
expressions  defining  the  prism  have  been  given.  These  enable  the 
manufacturing  specifications  of  a  prism  for  any  Bragg  cell  to  be  readily 
defined,  taking  full  account  of  its  Bragg  angle  and  polarisation  requirements. 
The  prism  would  also  be  useful  in  other  bulk  AO  and  Opto-SAW  architectures. 
The  effect  of  rotation  of  the  prism  has  been  well  defined.  It  causes  a 
variation  in  divergence  of  the  output  beams.  Because  of  this,  a  prism  made  to 
a  given  specification  can  be  used  with  a  range  of  different  Bragg  cells.  The 
variation  in  beam  divergence  can  be  calculated  from  the  analytic  expressions 
provided.  The  essential  requirement  that  the  crossover  point  of  the  input 
beams  lies  between  the  two  Bragg  cells  has  been  incorporated  into  the  design 
of  the  prism  and  the  relevant  equations  defining  crossover  distance  have  been 
provided. 

The  analysis  of  the  correlator  includes  expressions  for  the  light  beams  and 
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provides  a  level  of  detail  not  reported  explicitly  in  the  literature.  This 
gives  a  clear  understanding  of  the  processes  involved,  particularly  in  the 
case  of  phase  modulated  signals. 

Experimental  results  obtained  from  a  prototype  built  from  available  lead 
molybdate  cells  confirmed  the  usefulness  of  the  design.  Longer  transit  time 
cells  can  be  used.  This  would  increase  the  maximum  lag  between  unknown  and 
reference  signals.  The  use  of  50  ys  or  70  ys  tellurium  dioxide  cells  in  this 
architecture  would  result  in  a  practical,  compact  time  integrating  correlator 
suitable  for  use  with  a  variety  of  signals.  It  would  have  a  large  potential 
processing  gain,  currently  limited  by  the  state  of  the  art  CCD  linear 
photodiode  arrays  but  capable  of  readily  taking  advantage  of  advances  in  this 
area. 
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NOTATION 

nth  term  in  LRS 
amplitude 
bandwidth 
unit  vector 
crossover  distance 
energy  in  signal  s(t) 

Fourier  transform  of  t(t) 
frequency 
carrier  frequency 
unit  vector,  x  direction 
Bessel  function  of  order  i 
unit  vector,  z  direction 
refractive  index  of  light 
outward  normal  to  surface 
inward  normal  to  surface 
path  length 
rad ion 

position  vector 

reflectance  of  ith  surface 

cross  correlation  function  of  a  and  b 

real  part  of  expression  within  braces 

integer,  switch  number 

signal 

Hilbert  transform  of  s(t) 
time 

transmittance 
pulse  duration 
transit  time  of  cell 


Ui(x,t)U2(x,t) 


light  intensity  distribution 
correlator  output  light  intensity 
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v 

V 

W 

x 

z 

z(t) 

Z(w) 

a  .  . 

ij 

a(t) 

a.e.r 

A 

e(t) 

9  1 . 9  2  >  0  3 

0  . 

1 

9B 

Q 

X 

X 

o 

T 

V 
* 

4> 

u 


velocity  of  sound  in  cell 
signal  intensity 
plane  width 

spatial  coordinate  parallel  to  sound  propagation  in 
Bragg  cell 

spatial  coordinate  perpendicular  to  sound  propagation  in 
Bragg  cell 

analytic  signal  corresponding  to  s(t) 
the  Fourier  transform  of  z(t) 
angle  between  surface  i  and  surface  j 
phase  angle 
general  angles 

variation  of  parameter  following 
the  envelope  of  z(t) 
general  angles 
incidence  angle 
Bragg  angle 
wavelength  of  sound 
wavelength  of  light  in  medium 
wavelength  of  light  in  vacuo 
t  ime 

frequency  of  light 
beam  divergence 
general  angle 
angular  frequency  =  2irf 
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APPENDIX  I 

DETERMINATION  OF  PRISM  CHARACTERISTICS 


I . 1  Prism  angles 

The  prism  has  two  horizontal  faces  and  five  vertical  faces  (see  figure  4). 
Specification  of  any  four  angles  between  faces  determines  the  prism.  The 
two  input  beams  are  shown  figure  11.  The  requirement  that  the  two  outgoing 
beams  (one  reflected,  and  one  reflected  and  refracted)  should  be  divergent 
at  four  times  the  Bragg  angle  8g,  is  used  to  determine  the  prism  shape. 

Two  rays  are  considered.  ABCDE  (see  figure  11)  enters  surface  3  normally 
at  the  point  B.  This  intercepts  surface  2  at  C  making  an  angle  Jf  with  the 
internal  normal  n2  .  It  is  reflected  towards  surface  1  intercepting  it  at  D 

and  is  refracted  out  of  the  prism  to  E  where  it  is  required  to  make  an 
angle  of  29  with  the  reference  line  through  D  parallel  to  the  top  edge  of 
surface  3.  The  angle  the  emerging  ray  makes  with  the  normal  to  surface  1 
is  denoted  by  6-  The  second  ray  considered  is  FGH  which  intercepts 
surface  1  and  G  and  is  reflected  to  H  making  an  angle  a  with  the  normal  to 
surface  1  and  is  required  to  make  an  angle  20  with  the  reference  line 
through  G  parallel  to  the  top  edge  of  surface  3. 

The  two  rays^considered  are  initially  parallel  so  that  the  angle  FGH  is  2a, 
FGH  =  tt/2,  IGH  =  28g  and  FGH  =  FGI  +  IGH  so  that 

a  =  tt/4  +  8g  (I.  1) 

The  ang^le  FGJ  =  tt/2  -  a  and  FGJ  +  DGI  =  tt/2  so  that  DGI  =  PDK  =  a.  Also 
PDK  +  KDE  +  EDnt  =  it/2,  which  implies  that 


8  =  tt/2  -  a  -  28g 


Using  equation  (1.1)  this  becomes 


8  =  tt/4  -  30g 


(1.2) 


Ay  /  A 

The  angle  CDn^  is  denoted  by  0.  Noting  that  JIG  =  tt/2  -  a  and  summing  the 

angles  of  the  triangle  QCD  an  expression  involving  I  can  be  written,  namely 
2Jf  =  it  -(tt/2  -  0')  -  (tt/2  -  a)  so  that 


y  =  Me'  +  a)  (1.3) 


Denoting  MCN  by  ?  and  CDL  by  \p  it  is  clear  that  CDL  =  DCN  and 
C  +  4>  +  JT  =  tt/2.  Noting  that  LDn^  =  riiDK  implies  tp  -  tt/2  -  8'  -  a.  The 

A 

angle  al2  between  surfaces  1  and  2  can  now  be  found.  al2  =  MCN  or  a  -  C  so 
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Substituting  for  and  2 f  and  noting  that  sin8/sin8' 
index  of  the  glass,  the  expression  for  a12  becomes 


=  n,  the  refractive 


al2  =  a/2  -  i  sin  1  (1/n  sinB) 

Substituting  for  a  and  8, 


(1.5) 


a12  =  it/8  +  i  8g  -  i  sin  1  [1/n  sin  (tt/4  -  38g)]  (1.6) 

The  angle  a14  is  determined  by  the  fact  that  the  ray  DO,  the  component  of 

CD  internally  reflected  should  be  normal  to  surface  4  to  minimise  strong 
light  in  the  prism.  Since  n'DO  =  8'  it  follows  that  ol4  =  8'  or 

alu  =  sin  1  (1/n  sin8)  so  that  substituting  for  8 


alu  =  sin  1  [1/n  sin  (tt/4  -  38  ) 

D 


(1.7) 


The  internal  prism  angle  a13  between  surfaces  1  and  3  is  simply 
o  j  j  =  it  -  o  or 


a i  3  =  3  tt/4  -  0. 


(1.8) 


Surface  5  is  arbitrary  and  is  defined  to  be  normal  to  surface  3,  ie 
o3S  =  ir/2.  The  length  of  surface  2  is  such  that  the  refracted  ray  at  G 

intercepts  surface  4  not  surface  2.  The  angle  a32  is  equal  to  2f  so  that 

from  3,  and  substituting  for  a,  8'  and  8 


=  J  sin  1  [1/n  sin  (tt/4  - 


30g)]  +  (if/8  +  *  8g) 


(1.9) 


—  TT_a21-a41 


As  a  simple  check  the  following  relationship  should  be  true  (see  figure  4) 


CT  3  2  +  “3  1  + 

“14  + 

“4  2 

=  360° 

refractive  index  of  n 
made,  the  angles  are 

=  1.51509 
calculated 

of  the  BK7 
to  be: 

“32  = 

36° 

17' 

28" 

“3  l  = 

134° 

35' 

56" 

alk  = 

27° 

10' 

52" 

“1,2  = 

161° 

55' 

44" 
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and  the  sum  is  found  to  be  360°  O'  0M  as  required. 

1.2  Phase  equalisation 

There  is  a  phase  difference  between  the  two  col  linear  first  order  beams 
+  1^,  +lg  (see  figures  1(a)  and  2)  reaching  the  TIC  correlator  detector. 

This  is  due  to  their  different  optical  path  lengths  both  in  the  prism  and 
in  the  cells.  Moving  the  prism  parallel  to  surface  2  effectively  alters 
the  prism  thickness.  (Length  of  normal  from  surface  2  to  surface  1  at  the 
point  under  consideration)  and  proportionally  changes  the  path  length  of 
one  beam.  However,  the  movement  required  for  2ir  change  in  phase  is  only 
X/32  (approximately)  or  20  nm.  This  is  therefore  not  a  practical  way  of 
equalising  phase. 

A  phase  plate  of  dimension  1  x  1  cm  x  0 . 2  mm,  coated  for  improved 
transmission  at  X  =  632.8  nm  was  used.  This  was  placed  in  the  beam 
directed  towards  surface  1.  If  the  width  of  the  plate  is  W  and  the 
rotation  angle  is  9  the  optical  path  change  Ap  is  given  by 

Ap  =  n'W(l  -  cos  0)/cos  8 

where  n  is  the  refractive  index  of  the  plate.  For  example  a  2  mm  thick 
plate  requires  a  rotation  1°  10'  36"  to  change  the  phase  of  X  =  632.8  nm 
light  by  2tt. 

1.3  Prism  rotation  functions 

Rotation  of  the  prism  relative  to  the  incident  light  ( ie  increasing  the 
angle  of  incidence)  changes  the  directions  of  both  outgoing  beams.  They  do 
not  vary  by  the  same  amount  resulting  in  a  change  in  the  divergence  of  the 
two  beams . 

Consider  a  rotation  of  the  prism  through  a  small  angle  A8  in  an 
anticlockwise  direction.  The  new  angle  of  incidence  on  surface  1  for 

the  reflected  ray  (see  figure  12(a))  is  given  simply  by  =  a  +  A0  and  the 

angle  02  which  the  outgoing  reflected  ray  makes  with  the  reference  (top 

edge  of  surface  3)  is  given  by  2a^  -  u/2  or  substituting  for  and  then  a. 


0!  =  28d  +  2A0  (I . 10) 

D 

The  ray  entering  surface  3  is  now  not  normal  and  suffers  a  refraction. 
From  figure  13(b)  the  following  relations  can  readily  be  obtained. 


=  tt/2  -  a, 2  -  A0' 

=  it/2  -  <;,  =  a, j  +  A0  ’ 

a  3  2  +  +  tn  =  tt/2 


Substituting  for  2Tn, 


can  be  written  as 
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C2  =  it/ 2  -  2  a32  -  A8 ' 


Also 

=  C2  +  a  +  tt/2 

so  that 

02  =  2  o32  +  A0'  -  ir/4  -  0D 

b 

The  outgoing  ray  makes  an  angle  02  with  the  reference  line.  Since 
0  2  +  B  2  +  a  =  tt/2 

02  =  ir/4  -  0_  -  sin  1  (n  sin[2  a32  +  A0 '  -  rr/4  -  0  ] } 

b  b 

Substituting  for  a32  from  equation  (1.9) 

2  oc32  -  ir/4  -  0  =  sin  1  [  1/n  sin(n/4  -  30)] 

D  b 

so  that  02  is  given  by 

02  =  tt/4  -  0  -  sin  1  (n  sin(sin  1  [1/n  sin(u/4  -  38„)]  +  A0 '  ) }  (I. 11) 

b  b 

where 

A8 '  =  sin  1  [1/n  sin  A8 j 

Note  that  when  A0  =  0 ,  A8 '  will  also  be  zero  and  the  equations  (I.  10)  and 
(I.  11)  reduce  simply  to 


0i  =  02  =  28g 


as  required. 

The  divergence  0  between  the  two  outgoing  rays  is  just  the  sum  03  +  02 . 
For  small  A8 ,  0  is  given  approximately  by  a  simple  linear  relationship 


0  =  CqA0  +  C3 

For  the  manufactured  prism  C  =0.78  (±.02)r  and  C3  =  0.028  r(=40  ). 
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APPENDIX  II 

REVIEW  OF  PERFORMANCE  PARAMETERS 
II. 1  Processing  gain 

Since  the  AO,  TIC  is  essentially  a  matched  filter  the  well  known 
assessments  of  matched  filter  performances (ref . 14 , 15)  are  applicable.  The 
system  performance  criterion  is  processing  gain,  PG,  defined  to  be  the 
ratio  of  the  output  signal  to  noise  ratio,  SNRq  to  the  input  signal  to 

noise  ratio  SNR^.  By  def inition(ref . 14) . 

<;nr  =  Pea^  instantaneous  output  signal  power 
o  output  noise  power 


which  for  a  matched  filter  can  be  shown  to  be(ref.l4),  2E/kNo>  where  E  is 

the  received  signal  energy  and  No/2  is  the  noise  spectral  density  (in 
watts/cycle/second).  The  constant  k  is  near  unity  and  will  be  considered 
equal  to  1  here.  Consider  the  noise  to  be  white  and  define  the  noise 
bandwidth  B  of  the  matched  filter  as  the  bandwidth  of  a  rectangular  filter 
with  the  same  maximum  gain  and  the  same  output  noise  power  as  the  matched 
filter.  Then 


_  average  signal  power  at  the  filter  input  _  E/T 
i  noise  power  in  the  matched  filter  band  BNq 

where  T  is  the  effective  duration  of  the  signal.  Thu  processing  gain  is 
given  by  PG  =  SNRQ/SNRi  =  2  TB . 

A  number  of  specific  treatments  of  PG  in  TICs  exist(ref.l6  to  19,  and  23). 
It  is  relevant  to  review  that  of  reference  17.  This  deals  with  the  case  of 
additive  stationary,  zero  mean  Gaussian  noise.  The  output  signal  to  noise 
ratio  is  defined  in  terms  of  the  mean  and  variance  of  the  autocorrelation 
function. 


R(x) 


=  ±f 

TUt 


[s( 


I 


t  -  x/V)  +  n(t)]  s(t)dt 


The  mean  is  the  maximum  value 
variance  is  found  to  be 


P  of  the  autocorrelation  function  and  the 
s 


var  =  E 


If 


n(t  +  t)s(t)dt 


so  that 


nr* 


mean/var 
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When  the  ratio  of  the  bandwidth  of  the  noise  to  the  bandwidth  of  the  signal 
is  small  (ie  for  a  narrow  band  jammer)  the  conventional  PG  is  obtained. 
Reference  18  quotes  the  above  definition  of  PG  and  derives  expressions  in 
terms  of  SNR^  and  time  bandwidth  product  for  a  number  of  correlators. 

These  types  of  expression  are  useful  for  comparison  purposes. 

I I. 2  Dynamic  range 

In  practice  only  a  fraction  of  the  large  time  bandwidth  capability  of  a  TIC 
can  be  used.  The  system  is  limited  by  the  dynamic  range  of  the  output 
detector  array(ref . 1 ) . 

The  Bragg  cell  has  an  inherent  dynamic  range.  This  is  determined  by  the 
difference  between  signal  intensity  and  noise  sources.  The  most 
significant  noise  sources  are  the  third  order  intermodulation  components 
2f^  -  f .  formed  from  frequencies  f^  and  f ^  which  may  be  present  in  the 

cell.  In  the  Bragg  regime  this  intermodulation  intensity  is  given(ref . 20) 
by  |J3  (V)  | 2  -  C V / 2 ) 6  where  the  intensity  of  both  signals  are  equal  to  V, 

and  J3  is  a  third  order  Bessel  function.  A  graph  showing  the  limits  on 

Bragg  cells  inherent  dynamic  range  for  this  case  is  shown  in  figure  7. 
Some  of  the  analytic  results  were  verified  using  a  Bragg  cell  mathematical 
model (ref .10,11). 

The  lower  limit  is  due  to  optical  noise  due  mainly  to  diffraction  at  the 
aperture  restricting  the  beam.  Other  optical  noise  sources  are  discussed 
in  reference  22.  If  N  is  the  number  of  resolvable  frequencies  (bandwidth 
divided  by  resolution,  =  Bt)  then  the  optical  background  is (ref. 21). 


I  •  /I 

nun  o 


(  1/ttN)  z 


For  this  prototype,  B  =  30  MHz,  t  =  3  ps,  so  that  Bt  =  90  and  the  optical 

background  =  10  5 .  For  a  more  representative  N  of  say  1000  this  would  be 

=  10  7 .  The  limiting  dynamic  range  is  that  of  the  detector  array.  This 
point  was  discussed  in  Section  6  of  the  main  text. 

II. 3  Resolution 


The  divergence  a  of  the  first  order  from  the  incident  beam  is  given  by 


where  f  and  v  are  the  frequency  and  speed  respectively  of  the  sound  in  the 
cell.  Varying  the  frequency  will  vary  a  so  that,  Ax  =  XAfg/v. 

If  a  light  beam  passes  through  an  aperture  D  there  is  an  unavoidable 

diffraction  spread  a  .  =  XD.  Since 

r  min 


maximum  angular  diffraction 
angular  resolution 


frequency  range 
frequency  resolution 
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Ao 


a  . 
min 


then 


6f  = 


_Af 
=  6f 


v/D  . 


In  a  modulator  the  beamwidth  D  is  small.  In  a  deflector  it  is  usually  the 
maximum  cell  aperture.  In  the  latter  case  the  transit  time  T  =  D/v  so  that 
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Figure  2.  Bulk  AO  TTC  correlator  architecture 
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Figure  6.  Signal  and  noise  spectra 
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Figure  8.  Correlator  output  for  an  00k  signal 
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Figure  11.  Prism  geometry 
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Figure  13.  Crossover  distance 
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SUMMARY  OR  ABSTRACT: 

(if  this  is  security  classified,  the  announcement  of  this  report  will  be  similarly  classified) 


This  report  describes  the  design  and  testing  of  a  prototype 
correlator  using  lead  molybdate  Bragg  cells.  It  indicates  the 
suitability  of  this  design  for  cells  made  of  materials  with 
longer  transition  times,  such  as  tellurium  dioxide.  Compactness 
and  simplicity  are  achieved  by  the  use  of  a  specially  designed 
prism:  An  analysis  of  the  AO,  TIC  process  and  performance 
parameters  is  given.  .  : 
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